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Production mechanisms of charm hadrons 

in the string model 



In the hadroproduction of charm in the context of string fragmentation, the 
pull of a beam remnant at the other end of a string may give a charm hadron 
more energy than the perturbatively produced charm quark. The collapse 
of a low-mass string to a single hadron is the extreme case in this direc- 
tion, and gives rise to asymmetries between charm and anticharm hadron 
spectra. We study these phenomena, and develop models that describe the 
characteristics not only of the charm hadrons but also of the associated 
event . 
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Abstract 



Asymmetries in the production spectra of charm and anticharm hadrons (or generally, 
asymmetries between leading and non-leading particles) have been observed since long 
]T||. But it is only in recent years that the precision has improved so as to allow 

more detailed studies, especially in 7r~p events. Perturbative QCD calculations predict 
only very small asymmetries [0, |5|, so the origin of the observed asymmetries has to be 
sought in nonperturbative physics. 

Furthermore, the observed longitudinal momentum spectra of charm mesons are about 
as hard as or, in some cases, even harder than the perturbatively calculated charm quark 
spectra. This runs counter to naive expectations, e.g. based on experience from e + e~ 
annihilation experiments in the 10 GeV region, where the charm hadrons only keep about 
two thirds of the original charm quark energy 0. This number ought to be relevant for 
the hadroproduced charm events, since the Q 2 scales involved are comparable — by QCD 
scaling violations, the fraction slowly drops with energy. 

Several scenarios have also been proposed to explain the phenomenology @, |^, ITj ] • The 
two most frequently used ones are probably those of string fragmentation || and intrinsic 
charm ||10|| . In this letter we develop and study the string fragmentation approach. The 
key objective is to study the collapse of low-mass strings into a single hadron, which 
provides the main mechanism of flavour asymmetries in our approach. The outline of 
the letter is as follows. First we classify the different production channels for D mesons 
in hadron-hadron collisions, from perturbative and nonperturbative perspectives. Then 
we identify the origin of the charm asymmetry and compare the current model with 
data. The critical aspects of charm production are thereafter studied in more detail and 
different model variations are proposed to better understand the data. Further details 
and applications are intended to appear in a future paper. 

The production of a charm hadron can be subdivided into two steps: first the pro- 
duction of a cc pair, followed by the hadronization of these quarks. We will assume that 
the first step of this process is adequately described by standard perturbation theory 
and conventional parton distributions, i.e. without the inclusion of any intrinsic charm 
component in the proton wave function. The charm sea is thus assumed perturbatively 
calculable and peaked at small x values. Intrinsic charm and other modifications to stan- 
dard perturbative results may well exist at some level, but here we want to show that it is 
possible to understand existing data without invoking new mechanisms for the production 
stage. 

Also with this restriction, it is not possible to obtain unambiguous perturbative pre- 
dictions: results are sensitive to ill-determined parameters such as the charm mass. As 
a simplification, we will stay with lowest-order matrix elements for charm production, 
augmented by a parton-shower approximation to higher-order corrections. Given the 
uncertainties already noted, this appears adequate. 

The interesting phenomenology thus appears in the hadronization stage of our model. 
Here the partons of the hard interaction and of the beam remnants are connected by 



strings, representing the confining colour field [II|. Each string contains a colour triplet 
endpoint, a number (possibly zero) of intermediate gluons and a colour antitriplet end. 
The string topology can usually be derived from the colour flow of the hard process. For 
instance, consider the process uu —> cc in a 7r~p collision. The colour of the incoming u 
is inherited by the outgoing c, so it will form a colour- singlet together with the proton 
remnant, here represented by a colour antitriplet ud diquark. In total, the event will thus 
contain two strings, one c-ud and one c-d (Fig. |l]a). In gg — > cc a similar inspection 
shows that two distinct colour topologies are possible. Representing the proton remnant 
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Figure 1: Examples of different string configurations in a 7r~p collision: (a) uu — > cc has 
a unique colour flow; (b,c) gg — > cc with the two possible colour flows. 

by a u quark and a ud diquark (alternatively d plus uu), one possibility is to have three 
strings c-u, c-u and d-ud (Fig. |I]b), and the other is the three strings c-ud, c-d and u-u 
(Fig. |l|c). In addition to the valence uu annihilation and gg fusion production mechanisms, 
there are of course other possibilities, for example events involving sea quarks. These are 
included in the model but give small contributions to asymmetries and other observables 
at fixed-target energies, and therefore will not be discussed specifically in the following. 

In a process with two (or more) allowed colour flows, such as gg — > cc, the amplitude 
for either of them is calculable in perturbation theory, but with a nonvanishing interfer- 
ence term between the two. This term, corresponding to an indeterminate colour flow, 
fortunately is suppressed by a colour factor 1/Nq = 1/9. Its contribution to the total 
charm cross section therefore can be split between the two well-determined colour flows, 



e.g. according to the pole structure of the terms |[L2| , without much resulting ambiguity. 
Furthermore, we neglect the possibility of the cc pair forming a string, either by a per- 
turbative or nonperturbative colour (re) arrangement. Such mechanisms are likely to play 
a significant role for J/ijj production, e.g., but should be less relevant for the open charm 
production to be considered here. 

Once the string topology has been determined, the Lund string fragmentation model 



TT| can be applied to describe the nonperturbative hadronization. Assuming that the 
fragmentation mechanism is universal, i.e. process- independent, the good description of 
e + e~ annihilation data should carry over. The main difference between e + e~ and hadron- 
hadron events is that the latter contain beam remnants which are colour-connected with 
the hard-scattering partons. The structure of these remnants is not calculable from first 
principles, so this introduces some arbitrariness not constrained by e + e~ data. In the 
present model these aspects are parameterized in beam remnant distribution functions to 
be considered later. 

Depending on the invariant mass of a string, practical considerations lead to the need 
to distinguish three hadronization prescriptions: 

1. Normal string fragmentation. In e + e~ collisions the string system has a mass equal 
to the full CM energy, neglecting the not-too-frequent g — ► qq shower branchings 
which splits a string in two. This situation is ideal for an iterative fragmentation 
scheme, for which the assumption of a continuum of phase-space states is essential. 
The average multiplicity increases linearly with the string 'length' which means, 
neglecting gluon-emission effects, logarithmically with the string mass. In practice, 
this approach can be used for all strings above some cut-off mass of a few GeV. 

2. Cluster decay. In a hadron collision at fixed-target energies it frequently happens 
that a colour-singlet system contains two partons moving in the same general di- 
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Figure 2: Strings (dashed) in a 7r~p collision corresponding to the colour flows in Fig. |l]a, 
b and c respectively; (a) uu — > cc and (b,c) gg — > cc. If e.g. the colour-singlet system c-u 
in (b) has a small invariant mass we will call it a cluster and hadronize it by the procedure 
described in the text. 

rection, see Fig. |2|. This gives the system a small invariant mass, for which maybe 
only two-body final states are kinematically accessible. The continuum assumption 
above then is not valid, and the traditional iterative Lund scheme is not applicable. 
We call such a low-mass string a cluster, and treat it differently from above. When 
kinematically possible, a c-q cluster will decay into a charmed meson and a light 
meson by the production of a light quark-antiquark pair in the colour force field 
between the two cluster endpoints, with the new quark flavour selected according 
to the same rules as in normal string fragmentation. The decay kinematics is for 
the time being assumed to be isotropic in the rest frame of the cluster, which is the 
behaviour to be expected in the limit of vanishing phase space. 
3. Cluster collapse. This is the extreme case of the above situation, where the string 
mass is so small that the cluster cannot decay into two hadrons. It is then assumed 
to collapse directly into a hadron resonance, inheriting the flavour content of the 
string endpoints. The original continuum of string/cluster masses is replaced by a 
discrete set of hadron masses, mainly D and D*. By local duality arguments [[0| 
we assume that this does not change the total rate of charm production. This is 
related to the argument used in the e + e~ — > cc channel, that the cross section in 
the J/ip and ip' peaks is approximately equal to a purely perturbatively calculated 
cc production cross section restricted to the below-DD-threshold region. Similar 
relations have also been studied e.g. for r decay to hadrons |T3[, and there shown 



to be valid to good accuracy. In the current case, the presence of other strings in 
the event additionally allows soft-gluon exchanges to modify parton momenta as 
required to obtain correct hadron masses. Traditional factorization of short- and 
long-distance physics would then also protect the charm cross section. Local duality 
and factorization do not specify how to conserve the overall energy and momentum 
of an event, when a continuum of cd masses is to be replaced by a discrete D~ one, 
however. This will therefore be one of the key points to be studied below. 
Nature will not be so crude in its classification as we have been here, but hopefully our 
ansatz should be close enough to give a good first approximation. 

Charm cross sections are often presented as functions of Feynman-x, xp = Ph/Ph,max ~ 
2pl/£'cm, i-e. the longitudinal momentum fraction of the meson in the rest frame of the 
event |2], . The xp distribution for the different production channels in the standard 
Lund model is shown in Fig. || both for D~ and D + . We use Pythia 6.1 [15] with 
standard parameters to generate the plots. 
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(a) D production 



(b) D + production 




Figure 3: D T meson production in a 7r~p collision at a ir~ beam momentum of 500 GeV, 
using the default Lund Model, xp distribution (normalized per cc event) of (a) D~ and 
(b) D + for different production channels: (i) Cluster collapse, light quark from p end, (ii) 
Cluster collapse, light quark from 7r~ end, (iii) Cluster decay, light quark from p end, (iv) 
Cluster decay, light quark from 7r~ end and (v) String fragmentation, (c) The resulting 
asymmetry, Eq ([!]). Also shown is data from 0, ||. [|. 



In perturbative QCD the xp spectra of produced charm/anticharm quarks are identical 
to leading order, and the effects of higher orders are very small in this respect [[5], |(|. 
Therefore any asymmetry between charm and anticharm hadrons is a simple measure of 
nonperturbative effects. In 7r~p experiments the charm asymmetry is traditionally defined 

as 

x a(D-)-a(D+) 

A = A ^ = „(D-) + „ { D + y (1 > 

which measures the relative abundance of D~ over D + as a function of xp and p±. We 
will here mainly consider the variation with xp, integrated over p±, since A does not 
vary much with p± over the experimental range of xp values || |J. The asymmetry in 
Pythia is compared to experiment in 0, |3], |[, and the model shows general agreement 
with experiment, except in the region < xp < 0.4 where the model is predicting a larger 
asymmetry than found in the data. One collaboration 0] has tuned the parameters of 
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Pythia to obtain a good description of data also in this region, see below. 

In our model the main source of asymmetry for all xp > is the production of D~ 
(cd) mesons via cluster collapse involving a valence d from the pion beam, Fig. |3|a. This 
production channel is not open for D + (cd) since the ir~ does not contain a valence d. A 
number of model aspects influence the rate and character of cluster collapses, and in the 
following we will review the main ones. 

The charm mass enters in the perturbative matrix element as well as in the phase 
space, and so has a strong influence on the total charm cross section, but we have checked 
that the leading-order xp distribution is not affected much. It also, together with the 
light- quark masses, sets the threshold of the charm-cluster mass spectrum. The c mass 
to use need not even be the same in the two applications, cf. the familiar distinction 
between current algebra and constituent masses. A large combined mass of the charm 
and light quarks in a cluster dramatically reduces the importance of the cluster collapse 
mechanism. 

In the valence uu annihilation mechanism, the ir~ beam remnant is a d quark, so no 
internal structure need be specified. In the gg fusion process, however, the remnant is a 
ud system in a colour octet state, i.e. attached to two strings. A convenient approach is to 
imagine this system split into two separate u and d string endpoints. The beam remnant 
distribution function (BRDF) is introduced to describe how the (light-cone) momentum of 
the remnant is shared between the two, in fractions x an d 1 — Xi respectively. For an octet 
meson remnant the \ distribution is always implicitly symmetrized between the q and q, 
while for an octet baryon remnant one quark (picked at random among the three) takes 
the fraction \ and the remaining colour antitriplet diquark 1— %. To study the dependence 
of the asymmetry on the BRDF, we will consider two extreme cases and one intermediate. 
In one extreme one quark tends to take a small fraction of the available energy, much like 
the parton distributions in a hadron. In the other extreme, naive counting rules are used 
and the energy is, on the average, shared evenly between the quarks in the remnant. In 
Fig. [| the individual spectra as well as the asymmetry is shown for the full xp range. 

It is interesting to note the difference between the regions xp < and xp > 0. For 
xp > the D™ is a leading particle and the asymmetry is attributed to cluster collapses 
involving a d-quark from the pion beam. In the proton fragmentation region (xp < 0) the 
D~ is still leading and D + non- leading, so naively you would expect the asymmetry to be 
positive. However, when using an uneven sharing of energy in the proton beam remnant, 
diquark effects become prominent. As is seen in Fig. ^| the c is always connected to 
diquarks and the c is connected to quarks. Therefore D + mesons produced from a cluster 
(containing a diquark) via cluster decay will, on the average, be harder than D~ from 
cluster collapses, cf. Fig. |3|. We see that there are two competing effects for xp < 0, one 
that favours D + (diquark effect), and one that favours D~ (cluster collapse). The strength 
of the diquark effect depends strongly on the assumed energy sharing in the proton beam 
remnant. This is seen clearly in Fig. [|c where an even sharing makes the asymmetry 
change sign for xp < 0. If this region could be examined experimentally (p beam on p 
target) it would give us a hint at which distribution should be used. Possibly one could 
use data from other (related) experiments to reasonable choice of BRDF. 

The partons entering the hard interaction are traditionally taken to have a nonvanish- 
ing primordial k±, seen as a purely nonperturbative Fermi motion of partons inside the 
incoming hadrons, in addition to whatever is provided by perturbative gluon emission. 
Typical values should thus be 300-400 MeV. In many connections, also for charm in the 
current energy range 0, it has been noted that much higher values are required, at or 
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(a) D production (b) D + production 
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Figure 4: xp distribution of (a) D _ and (b) D + for different production channels (cf. 
Fig. [| which corresponds to the default of extreme uneven sharing) in a 7r~p collision at 
a 7r~ beam momentum of 500 GeV using an even sharing of energy between the quarks in 
a beam remnant, (c) the resulting asymmetry using three different choices of BRDF's: oc 

{ l ~x) k /\lx 2 + i, oc (l-x)7 \jx 2 + Cmin and oc (1-x)' respectively, with (k,l) = (1,0) 
for the pion and (3, 1) for the proton remnant, and c min = 0.6 GeV/E cm . They correspond 
to an uneven (full), intermediate (dashed) and even (dotted) sharing of energy. For the 
tt~ the three cases correspond to < x > = 0.14, 0.23 and 0.5 respectively. Also shown is 
data from [|, g, [§ . 

above 1 GeV. This remains somewhat of a mystery, which we do not attempt to solve 
here. The choice of primordial k± distribution is of non-negligible importance, both by 
providing a p± kick to the produced charm quarks and, by momentum conservation, an 
opposite kick to the beam remnants. 

The cluster mass spectrum is affected both by the choice of charm mass, of BRDF's 
and of primordial k± distributions. As a typical example, the c-d singlet produced in 
Fig. |2|a can have the mass distribution of Fig. |5|a when calculated on the parton-level with 
m c = 1.35 GeV and m^ « m u « 0, or as in Fig. [5|b if some parameters are varied. But we 
know that the observed mass spectrum of produced particles consists of peaks at the D~ 
and D*~ masses and then a continuum above the D~7r° threshold, cf. Fig. ||a. 
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Figure 5: (a) Distribution of cluster (full) and meson (dashed) masses in the string model. 
Clusters within the gray area collapse to D or D*~. (b) Dependence of the parton-level 
mass distribution on some parameters of the model. 

The onset of the continuum depends on the assumed threshold behaviour. In one 
extreme, a Dtt pair is always formed when the cluster mass allows it. In another, this 
continuum state is only reached after a succession of thresholds: D71", D*7r, Dp, D*p, etc. 
Thus a fraction of events may collapse to a single resonance also above threshold, cf. 
the gray area in Fig. |^a. In the default implementation, exactly one attempt is made 
to form two hadrons by a random choice of allowed flavours and spins (using the string 
fragmentation relative probabilities), and is allowed to succeed if the sum of hadron masses 
is below the cluster mass. By increasing the number of attempts made before giving up, 
the behaviour interpolates to the first extreme above. Thus we see that the mass spectrum 
of collapsing clusters is affected in the low end by quark masses and in the high end by the 
transition to two-body states, altogether giving a large range of possible cluster collapse 
rates. 

When a collapse occurs, confinement effects have to project the continuum of string 
masses onto the observed hadron mass spectrum. Because of the aforementioned local 
duality and factorization arguments, the total area of the spectrum should be conserved 
in the process. How the projection should be done is not known from first principles, 
however. 

One conceivable strategy could be to introduce a weight function consisting of S func- 
tion peaks at the D~ and D*~ masses, with suitably adjusted normalizations, and then 
a step function at the D~7r° threshold. This weight function, when multiplied with the 
partonic mass spectrum, should then give the hadron-level mass spectrum. Such an ap- 
proach is not well suited for Monte Carlo simulation, since the string mass is a complicated 
function of a number of variables and therefore the S function cannot easily be integrated 
out. However, on general grounds, we do not expect the overall distribution of event 
characteristics to differ significantly between events with a c-d string mass exactly equal 
to the D~ one, and events where the string mass is maybe 100 MeV off. An appealing 
shortcut therefore is to accept all partonic configurations and thereafter introduce some 
'minimal' adjustments to the kinematics to allow hadrons to be produced on the mass 
shell. Such a strategy would be consistent not only with local duality arguments, but 
also with the presence of soft final-state interactions, i.e. the exchange of nonperturbative 
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gluons that can carry some amount of momentum between the low-mass string and the 
surrounding hadronic system. In the following we will therefore adopt the language of 
'gluons' transferring energy and momentum between the strings in a collision, while leav- 
ing unanswered the question on the exact nature of those 'gluons'. Specifically, we will 
not address the possibility of changes in the colour structure of events by such 'gluons'. 

In this letter we will consider two different choices of energy shuffling schemes that 
can be said to be of opposite nature. We will further show that, as far as the asymmetry 
is concerned, the observable differences are small. 

The first approach is the standard one in Pythia. It consists of shuffling momentum 
to the parton i in the event that has the largest invariant mass when combined with the 
cluster cl that is supposed to collapse to a D (say). The new i and D momenta are given 
by 

p D = (l - S)pd + m , 

p\ = (1 - e) Pi + 5 Pcl , (2) 

where e and 5 are determined by the conditions p\, = and pf — p\ = mf. Using 
the parton furthest away from the cluster has the advantage of minimizing the required 
momentum transfer, but does not offer a particularly appealing picture physically. In the 
following we therefore try to formulate a scheme where energy and momentum is shuffled 
to partons or strings in the vicinity of the cluster. 

This problem can be approached in several ways, and we want to mention a few. The 
basic idea is to emit a soft gluon from the cluster, taking away the energy and momentum 
needed to put the remainder on the hadronic mass shell. One way to do this is simply 
to rescale the four-momentum of the cluster in order to give it the correct invariant 
mass. The gluon can then have both positive and negative masses and energies, which 
may be conceptually unappealing but not forbidden in principle. Worse is that such a 
rescaling can give xp > 1, i.e. kinematical inconsistencies. Another approach is to let the 
cluster decay isotropically into a hadron and a massless gluon. The gluon can still have 
both positive and negative energies, but this is not a problem in itself since the gluon is 
connected to the nearest string piece, which is then hadronized according to the standard 
string-with-gluons scheme [16| . A gluon with positive energy will increase the energy and 
mass of the string while a gluon with negative energy will decrease it. It is the string 
which is the physical entity, and it of course has to have a positive energy and mass. 
The only observable effect of the extra gluon is a slight increase (decrease) in the average 
multiplicity in the phase space neighbourhood of the inserted positive (negative) energy 
gluon, but this effect is too small to be seen in an actual experiment. Even so, there are 
problems with the simple approach. For instance, if a negative-energy gluon is connected 
to a string that already has a small mass, it can give the string a negative mass and this 
is not permitted. So, if there are no sufficiently massive strings left in the event, the gluon 
cannot be connected anywhere, and this method cannot be used. Only if the radiated 
gluon has positive energy can the mechanism always be used without problems. 

In practice, therefore, the negative-gluon approach is too error-prone. When the clus- 
ter mass is smaller than the intended hadron mass another approach is instead used to 
simulate the same effect, i.e. to take up energy from the vicinity of the cluster: a nearby 
string is allowed to emit a gluon which, when absorbed by the cluster, gives it the right 
mass. In more detail, imagine a string with endpoints q and q. Form a weighted sum of 
the endpoint momenta 

p s = ap q + f3p^= p q + Pq , (3) 

PqPq PqPq 



8 




Figure 6: Some results for different cluster collapse algorithms: (a) The momentum spec- 
trum for D~ mesons produced from clusters (normalized to the total D~). The two top 
curves are with default parameters and the lower ones are with an increased charm mass 
(from 1.35 to 1.7 GeV). (b) the resulting asymmetry. 

so that the end of the string that is closest to the cluster is weighted up relative to the 
one further away. Thereafter define 

Pd = Pel + 5p s , (4) 

with 5 determined by the constraint = m^. The meson will then have the correct 
mass and the string endpoint momenta are scaled down by factors 1 — 5a and 1 — 5(3, 
respectively. Special cases need to be introduced to avoid e.g. 1 — 5a < 0, but these affect 
only a small fraction of the events, and can be left aside here. 

In summary, the new algorithm transfers a 'gluon' from the cluster to the nearest 
string if the original cluster mass is above the D meson mass, and transfers a 'gluon' 
the other way if the mass instead is below. Some results are shown in Fig. |[ As can 
be seen, the largest difference is in the distribution of the collapsed mesons at high xp, 
where the spectrum is somewhat softer in the new scheme than in the old one. Another 
difference is around xp = 0, where the old algorithm tends to push the meson away from 
the middle, creating an unphysical dip in the distribution. This effect is seen more clearly 
in the distributions for a larger charm mass. The distribution with m c = 1.7 GeV is also 
included to show the relative importance of the collapse mechanism on the asymmetry. 
Clearly the charm mass — by regulating the amount of collapses — provides a much 
larger uncertainty (except possibly at very large xp) than the choice of collapse scheme 
but, if nothing else, the new algorithm remedies some cosmetic problems of the old one. 
The new algorithm also changes some other characteristics of the event, e.g. the mean 
multiplicity is changed by ~ 2 %, but this too is a small effect compared to the other 
uncertainties that we have mentioned. 

The E791 collaboration has produced a tuned version of Pythia The specific 
parameters involved have already been discussed above, so we here only provide a brief 
summary: 

• The charm quark mass is increased from 1.35 to 1.7 GeV. The fraction of cluster 
collapses is thereby reduced by about 35 %, which reduces the asymmetry consid- 
erably for xp > 0, as demonstrated in Fig. |6|. This large charm mass may appear 
unrealistic but it is also possible to reduce the number of cluster collapses by in- 
creasing the light-quark masses and/or increasing the probability for a cluster with 
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Figure 7: The D 1 * 1 asymmetry as a function of (a) xp and (b) p\ for the modified version 
of Pythia described in the text: new cluster collapse mechanism, constituent masses for 
the quarks (m c = 1.5 GeV, = m u = 0.33 GeV), intermediate BRDF's, primordial 
k± = 1.0 and using two attempts to form two particles from a cluster, (a) is for all p± 
and (b) is for —0.2 < xp < 0.8. Experimental data is from |2], |3|, ||. 

a mass above the D7r threshold to decay into two particles. By using constituent 
masses both for the light quarks and the charm quark (e.g. m c = 1.5 GeV and 
md = m u = 0.33 GeV) and by allowing two attempts to form two hadrons from 
a cluster, the number of collapses is reduced by a comparable amount. A charm 
mass of 1.5 GeV is also suggested from calculations of total cross sections for charm 
production in next-to-leading order QCD ||, and given the theoretical uncertainty 
regarding quark masses this value is not unreasonable. 

• A somewhat less peaked BRDF, like our intermediate scenario. The asymmetry is 
reduced somewhat in the region < xp < 0.4, Fig. [|. A photoproduction experiment 
[ P"7|j also seems to indicate flat BRDF's, but these are not determined by any basic 
principles, and are poorly known, so further studies will be needed here. 

• The width of the Gaussian primordial k± distribution is increased from 0.44 to 1.0 
GeV. As we have noted above, such a number is unexpectedly large, but in agreement 
with other data, and therefore rather standard these days. This allows cluster 
collapses between a charm quark and a beam remnant to occur also at fairly large 
p± values, thus leading to an essentially ^-independent asymmetry. Additionally, 
the p± kick added to charm quarks and beam remnants tend to increase the average 
invariant mass of the produced clusters, thereby reducing the number of cluster 
collapses. 

Taken together, these parameter changes gives a good fit to data, both for the asymmetry 
and for the shape of the individual charm meson spectra. Fig. [7] shows the asymmetry 
using default Pythia with the following modifications: new cluster collapse mechanism, 
constituent quark masses, intermediate BRDF's, an increased intrinsic k± and with two 
attempts to form two hadrons from a cluster. The amount of cluster collapse is here 
reduced by about a factor of three compared with the default, to ~ 16%. 

In summary, we have in this letter described the string fragmentation approach to 
charm production in hadronic collisions. A number of uncertainties have been identified 
and studied in detail, in particular the transition from a continuous string-mass distribu- 
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tion to a discrete hadron-mass one. The conclusion is that the model can describe much 
of the existing data, but also that these data do not fully constrain the choice of model 
parameters. Further data on charm production in 7r~p collisions may provide further 
information, as may charm production e.g. in ep collisions. Applications include, among 
other aspects, bottom production in hadron colliders where the asymmetry between B 
and B mesons has to be understood in order to facilitate a study of CP violation effects 



lq|. We intend to return to these and other related topics. 
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